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ABSTRACT 


The  thermal  boundary  resistance  between  a  YBa2Cu307  thin- 
film  and  an  MgO  substrate  was  measured  under  conditions  of 
varying  heat  flux.  Heat  flux  was  varied  in  a  manner  to  explore 
any  hysteresis  effects  present.  It  was  concluded  that 
hysteresis  effects  are  present  and  are  most  likely  attributed 
to  changes  in  the  peeling  or  compressive  stresses  in  the  thin- 
film.  The  changes  in  the  peeling  stresses  may  not  be  fully 
relieved  after  cycling  of  the  heat  flux  or  may  have  caused 
micros t rue tural  changes  near  the  interface  resulting  in 
changes  in  microscale  heat  transfer  characteristics. 
Additionally,  finite  difference  method  were  used  to  model  the 
physical  situation.  It  was  found  that  boundary  resistance 
values  generated  by  the  computer  program  were  several  orders 
of  magnitude  less  than  experimental  value.  It  was  concluded 
that  finer  meshes  must  be  used  in  order  to  increase  the 
accuracy  of  the  results.  It  was  recommended  that  the  modeling 
be  redone  on  a  main  frame  computer  using  finite  element 
methods . 
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I .  INTRODUCTION 


The  thermal  boundary  resistance  between  a  super¬ 
conducting  film  and  the  substrate  can  strongly  affect  the 
thermal  design  of  devices  made  from  high-temperature 
superconductors  (HTS) .  Heat  dissipation  has  always  been  an 
important  factor  in  the  design  of  electronic  devices.  In 
thin- film  HTS  devices,  poor  heat  transfer  may  result  in  the 
device  being  driven  from  the  superconducting  to  the  normal 
state.  In  a  superconductivity  application  this  would  prevent 
the  device  from  functioning  as  designed. 

The  current  use  of  superconducting  electronic  devices 
falls  into  two  basic  classes:  SQUID  systems,  which  are  used 
to  directly  measure  magnetic  flux;  and  Josephson  devices, 
which  perform  traditional  electronic  functions  [Ref.  1] .  The 
thermal  stability  of  thin  films  is  crucial  due  to  their 
critical  temperature  requirements  and  small  size,  and  is 
highly  dependent  upon  the  thermal  boundary  resistance  (Rb) 
between  the  film  and  the  substrate  [Ref.  2].  While  some 
characterization  of  Rb  has  been  done,  including  the  effects 
of  heat  flux,  further  investigation  is  warranted  [Ref.  2,3] . 
In  particular,  no  study  has  been  done  to  date  on  the 
hysteresis  effects  on  Rb  resulting  from  varying  the  heat 
flux. 
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Accordingly,  the  goal  of  this  study  is  to  vary  the  heat 
flux  applied  to  a  YBa2Cu307  thin  film  on  a  MgO  substrate  in 
the  normal  state  and  investigate  the  effects  on  Rb.  Heat 
flux  will  be  varied  at  several  different  substrate 
temperatures  and  will  be  cycled  to  observe  any  hysteresis 
effects . 

Additionally,  efforts  will  be  made  to  use  finite 
difference  numerical  techniques  to  model  the  physical 
situation.  Current  experimental  techniques  do  not  allow  the 
measurement  of  the  temperature  in  the  substrate  directly 
below  the  thin- film;  this  temperature  must  be  solved  for 
using  the  Laplace  equation.  Finite  difference  modeling  will 
yield  a  temperature  distribution  that  may  allow  refinement 
of  experimental  values  for  increased  accuracy. 
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II .  BACKGROUND 


A.  THERMAL  BOUNDARY  RESISTANCE 

1.  Heat  Transfer 

Whenever  there  is  a  contact  of  two  surfaces,  a  thermal 
resistance  to  heat  transfer  Rb  occurs  at  the  interface 
resulting  in  a  temperature  discontinuity.  This  is 
illustrated  in  Figures  1  and  2.  It  is  attributed  to 
imperfect  contact,  and  acoustic  and/or  lattice  mismatch. 

The  heat  transfer  across  an  interface  can  be  considered 
as  a  net  flow  of  phonons  (and/or  electrons)  between  the  two 
surfaces.  The  phonon  (and/or  electron)  transmission 
probabilities  can  be  calculated  to  determine  the  heat  flux. 
At  low  temperatures,  Rb  is  explained  by  considering  an 
acoustic  mismatch  between  the  two  solids  at  the  interface. 

If  the  acoustic  impedances  of  the  two  media  are  very 
different,  then  phonon  scattering  occurs. 

An  idealized  case  of  a  solid-helium  interface,  thermal 
boundary  resistance  is  often  called  the  Kapitza  resistance, 
which  is  explained  according  to  the  acoustic  mismatch  (AMM) 
theory  [Ref.  4].  This  theory  has  been  extended  to  solid- 
solid  boundaries,  however,  its  prediction,  except  at  low 
temperatures,  is  lower  than  empirical  values  [Ref.  5] . 
Recently,  the  AMM  theory  was  modified  to  include  the  effects 
of  anharmonic  interactions  in  crystals.  Inclusion  of 
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Figure  1.  Interface  Between  Two  Surfaces  [Ref.  6] 
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Figure  2.  Interface  Temperature  Discontinuity  [Ref.  7] 


these  Umklapp-processes  in  the  AMM  theory  (U-AMM)  resulted 
in  Rb  values  that  lie  closer  to  experimental  data  [Ref.  8] . 
However,  the  disparity  between  theoretical  and  experimental 
values  at  temperatures  >  30  K  still  necessitates 
experimental  characterization  of  Rb. 

2.  Characterization  of  Rb 

The  thermal  boundary  resistance  (Rb)  between  thin- film 
HTS  and  their  substrates  has  been  measured.  The  experimental 
measurement  techniques  are  fully  described  in  Phelan  et  al . 
[Ref.  9]  and  are  based  on  the  methods  introduced  by  Swartz 
and  Pohl  [Ref.  10] .  Using  this  technique,  previously 
measured  values  for  Rb  for  an  Er-Ba-Cu-0  thin-film  on  MgO 
substrate  with  a  Tc  of  approximately  69  K  were  found  to  be 
on  the  order  of  10'3  K-cm2/W  [Ref.  11]  . 

Previous  measurements  were  conducted  using  two 
different  etching  patterns  for  the  film:  two  parallel  strips 
of  equal  length  in  a  straight  pattern,  and  two  strips  of 
unequal  length  in  a  meander  pattern  (Figures  3  and  4) .  The 
straight  pattern  has  the  advantage  of  being  simple  and, 
thus,  easier  to  etch.  The  meander  pattern  is  more 
complicated,  but  was  designed  with  a  smaller  heater/sensor 
separation  than  the  straight  pattern,  and  was  meant  to 
provide  more  film  area  for  radiative  property  measurements. 
Both  patterns  have  only  one  heater  and  one  sensor  strip 
that,  if  damaged  during  fabrication,  will  render  the  sample 
useless . 
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Figure  3.  Straight  Pattern 


Figure  4 .  Meander  Pattern 


The  film  under  test  in  this  study  utilized  four 
parallel  strips  etched  in  the  YBa2Cu307  thin- film,  as  shown 
in  Figure  5 .  Should  an  etching  imperfection  occur  in  one  of 
the  strips,  other  adjacent  pairs  will  be  available  for  use. 
Both  strips  are  of  equal  width  and  either  can  function  as 
the  heater  or  the  sensor.  The  question  remains  whether  two 
non-adjacent  film  strips  can  be  used  for  data  acquisition. 
The  use  of  the  finite  difference  numerical  technique  may 
hold  the  answer. 

3 .  Mathematical  Modeling 

The  finite  difference  technique  will  be  used  to  model 
the  film- substrate  response  to  a  temperature  or  a  heat  flux 
imposed  on  one  of  the  four  thin- film  strips.  The  control 
volume  approach  will  be  used  to  obtain  a  graded  network  that 
can  be  varied  to  minimize  the  overall  number  of  nodes  in  the 
network,  yet  maximize  accuracy  within  the  heat  affected  zone 
[Ref.  12] .  Due  to  the  large  number  of  nodes  expected,  a 
Gauss-Seidel  algorithm  will  be  used  for  the  program  [Ref. 

13]  . 

The  results  of  the  program  will,  at  the  very  least, 
yield  a  qualitative  representation  of  the  temperature 
distribution  within  the  substrate.  The  goal  is  to  obtain  a 
value  for  the  temperature  of  the  substrate  directly  beneath 
the  heater  strip.  This  temperature  will  be  used  to  bias  the 
experimental  results  to  obtain  a  more  accurate  value  for  Rb. 
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B.  YBCO  THIN- FILM  SUPERCONDUCTORS 

1 .  YBCO 

Sintered  samples  of  YBa2Cu307.x  exhibit  superconductivity 
at  temperatures  above  the  boiling  point  of  liquid  nitrogen 
(77  K) .  The  structures  of  the  YBa2Cu307  unit  cell  and 
multicellular  YBa2Cu307  are  shown  in  Figures  6  and  7 
respectively.  The  Cu-0  planes  are  assumed  to  be  the  focus  of 
the  superconductivity  mechanism.  The  oxygen  concentration 
determines  whether  the  material  is  superconducting  or  not. 
The  onset  of  superconductivity  occurs  at  a  value  of  x=0.55 
(7-x=6 . 45)  as  the  YBCO  structure  changes  from  tetragonal 
(non  superconducting)  to  orthorhombic  (superconducting) . 

[Ref.  17] 

2 .  Thermal  Properties 

The  thermal  conductivity  (k)  of  YBa2Cu307  figures 
prominently  in  both  the  error  analysis  and  the  finite 
difference  methods  used  in  this  study.  Given  the  structure 
of  the  unit  cell,  it  is  not  surprising  that  the  thermal 
conductivity  is  anisotropic. 

Values  of  temperature  vs  thermal  conductivity  values  of 
YBCO  are  shown  in  Figure  8  [Ref.  18].  Thermal  conductivity 
in  the  c-plane  is  critical  to  the  heat  transfer  from  the 
heater  strip  into  the  substrate.  Thermal  conductivity  in  the 
ab-plane  figures  prominently  in  the  error  analysis  when 
evaluating  for  conduction  losses  in  the  film  lead. 


9 


o  o  o  o 

CM  rf  <D  00 


o 

o 

o 

o 

o 

o 

O 

O 

O 

O 

o 

o 

CM 

CD 

CD 

o 

CM 

CO 

CD 

O 

*■* 

T— 

T*“ 

CM 

CM 

CM 

CM 

CM 

CO 

T  (K) 


Figure  8.  YBCO  Thermal  Conductivity 


11 


Values  for  temperature  vs  thermal  conductivity  for  MgO 
are  plotted  in  Figure  9  [Ref.  19] .  Thermal  conductivity  for 
the  MgO  substrate  is  important  not  only  for  the  error 
analysis,  but  also  in  the  finite  difference  programs. 
Previous  studies  failed  to  use  accurate  values  for  thermal 
conductivity  for  both  YBCO  and  MgO  in  their  error  analysis. 
This  study  utilizes  cubic  spline  interpolation  to  the 
maximum  extent  possible  to  yield  the  best  values  for  the 
thermal  conductivities. 
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III.  EXPERIMENTAL 


A.  MATERIALS 

The  sample  used  for  this  study  was  an  epitaxial 
YBa2Cu307  thin  film  on  a  single-crystal  MgO  substrate.  There 
are  three  steps  in  the  fabrication  process:  film  deposition, 
microfabrication  by  photolithography,  and  gold  pad 
deposition  for  lead  connection. 

The  films  used  in  this  study  were  purchased  from  Excell 
Superconductor.  The  YBa2Cu307  film  was  deposited  on  the 
substrate  by  rf -magnetron  sputtering.  This  sputtering 
process  is  fully  described  in  Adachi  et  al .  [Ref.  20]  The 
MgO  substrate  thickness  was  1  mm  and  the  area  was  10  mm  x  10 
mm.  The  YBa2Cu307  film  was  3000  A  (300  nm)  thick.  The  film 
exhibited  a  critical  temperature  of  Tc  -  78  K  during 
testing. 

Etching  was  done  by  photolithography  using  a  pattern 
drafted  in  AutoCAD  10  by  polylines  (Figure  10) .  The  central 
region  in  Figure  10  was  a  result  of  the  close  spacing 
between  the  heater  and  sensor  strips  which  caused  them  to 
appear  to  merge  into  one  dark  band.  The  AutoCAD  pattern  is 
converted  into  a  mask  by  a  pattern  generator.  The  mask  is  a 
glass  plate  with  the  coated  circuit,  which  resists  the 
passage  of  ultraviolet  light.  Photoresist  is  spun  onto  the 
films  and  the  mask  steps  the  required  circuit  pattern 
optically  onto  the  film  [Ref.  21].  A  standard  developer  is 
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then  used  to  develop  the  photoresist  to  reveal  the  areas  of 
the  film  to  be  etched.  The  etching  solution  is  1%  H3P04 
(hypophosphoric  acid)  and  the  film  is  soaked  for  15  minutes 
without  agitation.  The  gap  between  the  heater  and  the  sensor 
strips  is  0.01  mm  which  made  etching  difficult,  however, 
high  quality  etching  was  accomplished.  Detailed  methods  for 
etching  YBCO  films  are  given  by  Shin  and  Qiu  [Ref.  22] . 

The  connection  pads  to  the  etched  film  were  deposited 
through  a  mask.  The  mask  was  fabricated  out  of  polished  and 
buffed  sheet  metal  brass  with  holes  to  match  the  etched 
pattern.  It  was  closely  pressed  onto  the  etched  film, 
adjusted  to  match  the  pattern  and  taped  in  place.  Gold, 

1000  A  thick,  was  then  deposited  to  form  the  connection 
pads.  The  resulting  etched  pattern  with  connection  pads  is 
shown  in  Figure  11.  A  16 -pin  connector  was  then  fabricated 
to  allow  electrical  connection  of  the  pads  to  the  data 
acquisition  unit  while  the  sample  was  in  the  cryo-cooling 
apparatus . 
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B.  EQUIPMENT 


To  prevent  contaminants  and  moisture  from  condensing  on 
the  sample,  experiments  are  conducted  at  a  very  high  vacuum 
(10'6  torr) .  Operating  at  a  high  vacuum  also  reduces  the 
thermal  mass  which  the  cryo-cooling  equipment  needs  to  cool 
and,  most  importantly,  drastically  reduces  the  heat  leakage 
from  the  sample  by  convection  in  the  chamber.  The  sample 
holder  containing  the  sample  is  placed  within  a  stainless 
steel  hood  (vacuum  chamber)  which  is  evacuated. 

Additionally,  a  brass  shroud  is  provided  between  the  sample 
holder  and  the  stainless  steel  hood  to  reduce  radiation  loss 
during  the  process .  The  seal  between  the  hood  and  the 
apparatus  is  accomplished  with  an  O-ring  and  Dow  Corning 
high  vacuum  grease.  Figure  12  shows  the  experimental  setup. 

A  CTI  Cryogenics  model  8001  (350  CP)  cryo-cooler 
is  used  which  can  produce  temperatures  down  to  19  K  under 
good  vacuum  conditions  (low  cooling  load) .  A  diffusion  pump 
is  used  in  conjunction  with  the  cryo-cooler  and  can  achieve 
vacuums  of  up  to  10  ^  torr  in  about  2-3  hours.  The  pump  is 
connected  to  a  roughing  pump  for  initial  evacuation.  The 
pump  is  cooled  by  water  and  supplemented  by  hand-poured 
liquid  nitrogen.  The  vacuum  level  is  measured  by  a  hot-wire 
•cathode  gauge.  The  diffusion  pump  set  up  is  shown  in  Figure 
13  . 
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Figure  12 .  Experimental  Setup 
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Figure  13 .  Diffusion  Pump  [Ref.  23] 


A  Keithley  current  source  was  used  for  the  heater  strip 
and  the  voltages  are  measured  by  a  Keithley  Nanovoltmeter . 
Previous  experimental  setups  utilized  a  two-phase  lock-in 
amplifier  from  EG&G  to  supply  a  small  ac  current  for  the 
sensor  strip  at  a  frequency  of  47  HZ  and  lock  onto  the 
signal  to  measure  voltage  drop  across  the  strip.  This  was 
found  to  introduce  noise  into  the  system  and,  consequently, 
the  lock-in  amplifier  was  removed  from  the  circuit.  All  dc 
voltages  across  shunt  resistors  were  measured  with  Hewlett 
Packard  voltmeters. 

In  previous  studies,  data  acquisition  was  accomplished 
using  a  Hewlett  Packard  computer  with  an  HP  basic  program, 
however, this  did  not  provide  the  operator  with  the  ability 
to  monitor  important  experimental  parameters  during  the  data 
taking.  Monitoring  and  data  acquisition  for  this  study  was 
accomplished  on  a  Power  Machintosh  8100/100  computer  linked 
through  an  IEEE-488  interface.  A  Labview  program  was  used 
to  provide  the  operator  with  a  visual  display  of  heater 
resistance  (ohms) ,  sensor  resistance  (ohms)  and  coldfinger 
temperature  (K) .  Additionally,  digital  readouts  of  heater 
current  (ma) ,  sensor  current  (ma) ,  heater  strip  voltage  (V) 
and  sensor  strip  voltage  (V)  were  displayed.  The  Labview 
program  enabled  better  control  of  critical  experimental 
parameters,  such  as  coldfinger  temperature,  than  were 
previously  achievable  and  allowed  the  operator  to  view  such 
events  as  the  onset  of  superconductivity  in  the  sample  at  Tc. 
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The  thin  film  deposited  on  its  substrate  is  mounted  in 
a  copper  plug.  The  copper  plug  acts  as  a  cold  finger  to 
conduct  heat  out  from  the  bottom  of  the  substrate.  The 
design  of  the  copper  cold  finger  is  shown  in  Figure  14. 
Apeizon-N  thermal  grease  is  used  between  the  substrate  and 
the  copper  block  to  ensure  good  thermal  contact .  The  pin 
connector  device  is  then  screwed  in  place  to  make  electrical 
connections  from  the  sample  to  the  data  acquisition 
equipment.  The  copper  cold  finger  plug  is  then  screwed  into 
the  metallic  chamber  which  is  connected  to  the  cryo-cooler. 

C .  PROCEDURE 

The  following  procedure  for  measuring  the  thermal 
boundary  resistance  (Rb)  between  the  film  and  the  substrate 
is  as  described  in  Phelan  et  al .  [Ref.  25] .  Referring  to 
Figure  15,  two  adjacent  YBCO  thin- film  strips  are  used:  one 
is  used  as  the  heater  and  one  is  used  as  the  sensor.  Four 
electrical  contacts,  including  the  voltage  taps  V*,.  and  Vf, 
are  made  to  each  strip.  This  allows  for  independent  four- 
point  measurement  of  the  resistance  of  each  strip  for 
various  bias  currents. 
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To  measure  the  thermal  boundary  resistance  between  the 
heater  strip  and  the  substrate,  a  relatively  large  dc 
current  is  passed  through  the  heater  strip,  while  a  much 
smaller  sensing  current  (Is5  juA)  is  passed  through  the 
sensor  strip.  The  Joule  heating  in  the  heater  strip 
increases  its  temperature  and  the  adjacent  substrate 
temperature,  and  thus  causes  the  temperature  of  the  sensor 
strip  to  increase.  Under  steady- state  conditions,  Rb  can  be 
easily  determined  from  the  definition  of  the  thermal 
boundary  resistance 


Rb= 


Th~Ts 


(1) 


where  Th  is  the  temperature  of  the  heater  strip  and  Ts  the 
temperature  of  the  substrate  directly  below  the  sensor 
strip.  Heat  flux,  q”,  is  accurately  calculated  from  the 
known  heater  voltage  and  current  along  with  its  contact  area 
(0.0016  cm2)  to  give  values  in  terms  of  W/cm2.  The 
temperature  of  the  substrate  below  the  heater  strip  is 
nearly  equal  to  the  temperature  of  the  sensor  strip,  since 
the  two  strips  are  only  separated  by  0.01  mm. 
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Both  Th  and  Ts  are  determined  from  the  steady- state 
values  of  their  electrical  resistances,  which  are  considered 
to  attain  steady  state  when  the  subsequent  data  points 
change  by  less  than  0.1%,  after  which  the  subsequent  10 
points  are  averaged.  After  all  data  acquisition  runs  are 
completed,  the  sample  is  cooled  to  well  below  Tc,  and  the 
cryo-cooler  is  turned  off.  Temperature  calibration  curves 
for  electrical  resistance  of  the  heater  and  sensor  strips 
are  determined  by  acquiring  data  as  the  sample  slowly  heats 
up  to  room  temperature.  Cubic  spline  interpolation  is  used 
to  extract  values  of  Th  and  Ts  from  the  temperature  and 
resistance  values  acquired  from  the  calibration  curve. 
Calibration  curves  for  one  set  of  data  are  shown  in  Figure 
16  . 
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IV.  RESULTS  AMD  DISCUSSION 


A.  Rb  CORRECTION 

Values  of  Rt,  are  calculated  and  tabulated  for  each  data 
run.  These  values  are  calculated  based  on  the  assumption 
that  the  temperature  of  the  substrate  below  the  heater  strip 
is  nearly  equal  to  the  temperature  of  the  sensor  strip.  The 
temperature  difference  for  each  data  run  must  be  corrected 
using  the  solution  to  the  Laplace  equation  [Ref.  26]. 

The  representative  heat  conduction  problem  is  shown  in 
Figure  17  and  the  solution  to  the  Laplace  equation  is  fully 
outlined  in  the  Appendix  I.  The  temperature  distribution  in 
the  substrate  is  given  by 

00 

T{x,y)  Cncos(3nx  sinh(3ny  (2) 

n-1 


where 


2gsinPnx1 

^p„acosh(3ny 


(3) 


and 


T(x,b) 


00 

-E 

n  =  l 


2gsin3nx1 

wU 


cos(3nx 


(4) 
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The  temperature  difference  in  the  substrate  between  the 
heater  and  the  sensor  strip  is  given  by 

Ar  A  2gsinp„x1  (  sinp^  sinp„x2 ^  (5) 

h  k$2na  xi  a_x2 


which  for  the  dimensions  of  the  sample  becomes 

AT=0. 015-2 
k 


The  temperature  in  the  substrate  directly  beneath  the  heater 
strip  is  found  by  adding  the  AT  calculated  from  equation  (6) 
to  the  measured  sensor  strip  temperature  (Ts)  : 

T  ~t  +Ar  (7) 

1  substrate  ±  s  ** 


Finally/  the  corrected  value  of  Rb  is  found  using 


^ substrate 


(8) 


Alternately,  Swartz  presents  another  method  for 
determining  Tsubstrate  [Ref.  27].  The  temperature  difference  in 
the  substrate  can  be  found  by  integrating  the  Laplace 
equation  using  Laplace  transforms. 
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The  result  is 


A  T= 


qd  r  (s-d)  ■ 


nk  2d2 


In  (s-d)  +  -fg+<*)  2  In (s+d)  -In (d)  -  — ln(s)  ] 
2d2  d2 


(9) 


where  d  is  the  width  of  the  heater  and  sensor  strips,  and  s 
is  the  separation  between  the  strips  from  center  to  center. 
The  resulting  AT  yields  corrected  Rb  results  that  are 
consistent  with  those  based  on  the  AT  calculated  from 
equation  (6)  . 

Applying  the  AT  obtained  from  equation  (6),  the 
corrected  Rb  is  calculated.  Figure  18  is  a  plot  of  Ts  vs 
experimental  and  uncorrected  values  of  Rb.  As  expected,  the 
corrections  are  small  and  yield  values  of  Rb  that  are 
slightly  less  than  the  experimental  values.  Uncertainty 
analysis  must  now  be  done  on  Rb  for  each  data  run. 

B.  UNCERTAINTY  ANALYSIS 

1.  Uncertainty  Equation 

Consider  Y  which  is  a  function  of  various  parameters 
represented  as  follows 

Y=f(xx,x 2,...,xn)  (10) 

where  x1 ,  x2, ...,xn  are  independent  variables. 


30 


The  uncertainty,  uf/  in  Y  can  be  written  as  [Ref.  28] 


Ur=  [  (  Uy 


AL)Z +u  (-M-)2+, 
dxx  ’  x*  dx2 


+Uy  ( 


df 

dx„ 


210.5 


)2] 


(ID 


where  ux1,  ux2,  .  .  . ,  uxn  are  the  uncertainties  in  measuring  x1 , 
x2,  .  .  . ,  xn,  respectively. 

The  uncertainty  for  Rb  is  given  by 


ur=1(u 


Ar' 


df 

dAT 


)  2  +  U  n ( 


df 
dg  ' 


)2] 


0.5 


(12) 


<?  g 

2 .  Temperature  Difference  Uncertainty 

The  temperature  difference  is  given  by  AT=Th-Ts,  and 
its  uncertainty  is  given  by 

uat=  [  (  utJ  2+  ( ur,)  2]  0-5  (13) 

The  temperatures  Th  and  Ts  are  obtained  from  the  calibration 
curves  for  the  electrical  resistances  of  the  heater  and  the 
sensor  strips.  Previous  studies  divided  the  calibration 
curves  up  into  sections  and  performed  a  polynomial  fit  on 
them  to  obtain  equations  for  the  curves.  The  uncertainty  was 
then  estimated  from  the  standard  deviation  of  the  difference 
between  the  calibration  data  points  and  the  fitted  curves. 
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The  present  study  uses  cubic  spline  interpolation  to 
determine  the  desired  resistance  value  from  the  calibration 
data,  thus  eliminating  the  need  for  curve  fitting.  The 
uncertainties  for  Th  and  Ts  are  therefore  reduced  to  the 
uncertainty  of  the  Si  diode  temperature  sensor  (±  0.1K).  The 
resulting  uncertainty  for  AT  from  equation  (13)  is  ±  0.14K. 

3 .  Heat  Flux  Uncertainty 

Uncertainty  in  the  heat  flux  results  from  uncertainty 
in  the  heat  transferred  from  the  heater  strip  into  the 
substrate.  The  heat  transferred  into  the  substrate  is  the 
measured  heat  transfer  minus  conduction  and  radiation 
losses : 

<^heater~(^measuied  ^cond  losses  ^rad  losses  (14) 


Its  uncertainty  is 


u  =[(u  )2+(u„  )  2  +  ( 11  -  )2]0-5 

Qbaatai  Qma&auxtKl  Qcond  Iosssm  Grad  losaaa 


(is: 


The  uncertainty  for  q^^^  is  equal  to  the  uncertainty 
in  the  power  applied  to  the  heater  strip  (q=EheaterIheater)  , 
which  is 


U^sur^[iU^I)2+{UIE)2] 


0.5 


(16) 


The  uncertainties  uE  and  u,  are  taken  as  the  uncertainties 
for  measuring  voltage  and  current  which  are  ±  1% . 
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The  uncertainty  for  qconduction  loss  is  expressed  as  the 
percentage  of  ^  which  is  conducted  through  the  film 

lead  to  the  substrate,  and  along  the  film  lead  and  through 
the  pin  connector.  Figure  19  shows  the  possible  heat 
conduction  paths  and  Figure  20  shows  the  thermal  circuit 
with  the  resistances  along  the  conduction  paths. 

Referring  to  Figure  20,  the  thermal  resistance  across 
the  YBCO  heater  film  is  defined  as  RfU(n  htr  and  the  resistance 
across  the  substrate  below  the  film  as  Rsubhtr-  The 
resistance  across  the  substrate  below  the  film  lead  is 
designated  Rsub  lead.  The  resistance  along  the  film  lead  is 
defined  as  RfUm  lead. 

The  conductive  thermal  resistances  are  given  by 


AT_  A T  _  A z 
q  q  ".A  k-A 


(17) 


where 

Az  is  the  thickness 
k  is  the  thermal  conductivity 
A  is  the  heat  transfer  area 
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Figure  19.  Heat  Conduction  Paths 
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Figure  20.  Thermal  Circuit 
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The  physical  dimensions  and  thermal  conductivities  used 
to  calculate  the  thermal  resistances  are  given  in  Table  1. 
The  ratio  of  the  heat  flows  is  determined  from  the  ratio  of 
the  total  resistances  in  the  three  branches  of  the  thermal 
circuit.  The  values  of  thermal  conductivities  and  branch 
resistances  are  from  analysis  of  a  typical  data  run  at  Ts  = 
91  K. 

The  value  for  RfU(n  lead  is  very  high  in  comparison  with 
the  other  resistance  values  and,  consequently,  the  dominant 
heat  leakage  in  the  film  lead  is  through  the  substrate. 

Thus,  referring  to  Figure  21,  we  need  to  determine  xcrit.cal  at 
which  the  ratio  qfUm  lead/substrate  and  qlead/pin  connector  is 
sufficiently  high  (taken  as  0.5). 

For  the  film  lead  of  length  x  and  width  w,  resistance 
between  the  film  lead  and  substrate  and,  the  resistance 
along  the  film  can  be  written  as 

J?  —  ^  (18) 

^ film  lead/ subtrate  t.  urir  v  ' 

KMgOwx 


R 


film  lead  ^^d 


(19) 
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Parameter 

Values 

Resistances 

Magnitudes 

^ybco  c-direction 

2.1496  W/m-K 

^  film.htr 

87.431  K/W 

^ybco  ab-plane 

9.0496  W/m-K 

Rb  x  area  under  film 

5.608  xIO^K/W 

^MgO 

295.66  W/m-K 

R  sub.htr 

2212  K/W 

Film  Thickness 

R  sub,  lead 

9915  K/W 

Length  of  Heater  Film 

Area  of  Substrate  Under  Heater 

R  film  lead 

19830  K/W 

Length  of  Film  Lead  (avg) 

^  pin 

0.002  K/W 

Width  of  Film  Lead  (avg) 

I 

Total  Area  of  Film  Lead 

Thickness  of  Substrate 

10  m 

Table  1.  Error  Analysis  Values 


Film 


Film  Lead 


^lead 

— 1 

1 - ri - 

i 

d 

I 


Figure  21.  Film  Lead  Heat  Leakage  Paths 
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This  yields  a  ratio  of  (,ead/substrate  ^ lead/pin  connector 


Qfilm  lead/substrate  _  X  ^ MgO 
Qfilm  lead 


kyBCoJ^t 


(20) 


Upon  assigning  a  value  of  0.50  to  the  ratio,  the  value 
for  x  .  becomes  2.143  mm,  which  gives  us  a  value  for 

critical 

Rsub,lead  of  19830  K/W  and  RfUn,  lead  °f  9915  K/W-  R1 '  R2  and  R3 
are  defined  as  the  total  branch  resistances  for  qheater/ 

<Jfi  Lm  lead/substrate  and  <3lead  respectively. 

The  ratio  of  the  heat  leakage  to  the  heat  supplied  to 

the  heater  is 

P  _  Qfilm  lead/  substrate*  Q lead/  pin  connector  _ _ ^1 _ _q  29  (21) 

theater  (  _i_  +  -1.  )  _1 

i?2  -^-3 


which  gives  a  percentage  heat  loss  due  to  conduction  of 


Qfilm  lead/  substrate*  Q lead/ pin  connector^  *100  =  —  *  39 

Qfilm  lead/  substate*  Qlead/  pin  connector* Qheater  1  •  ^ 9 


•100  =  22.7%  (22) 


It  is  reemphasized  that  the  values  for  Rb  and  the  thermal 
conductivities  used  in  this  calculation  were  from  a  typical 
data  run.  Individual  data  points  will  be  evaluated  using 
their  specific  information  to  obtain  their  uncertainty 
values . 
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4.  Radiation  Loss  Uncertainty 

Part  of  the  heat  from  the  film  radiates  to  its 
surroundings,  which  is  the  cold  head  holder.  The  cold  head 
holder  is  assumed  to  be  at  Tdiode  and  the  resistance  due  to 
radiation  heat  transfer  is  determined.  This  evaluation  is 
done  using  the  Stef an-Boltzmann  law  to  calculate  the  heat 

flux  between  ^heater  anC*  ^diode* 

g//=e  a  [  T^eater~Tdiode^  ^ 2  3 


where 

q"  =  heat  flux  from  the  film  to  the  cold  head  holder 
6  s  emissivity  of  the  film  taken  as  1 
a  =  Stef  an-Boltzmann  constant  (5.6697  x  10  8  W/m2-K4) 
Expanding  the  term  [T4heater-T4diode]  in  equation  (23), 

qr//=eo  [T2hCl-T2diode]  [T  ACr+T  diode ^ 
q"-tO  [ThtI+Tdiodg]  [ ThtI-Tdiode ]  [T  hti  +  T  diode ^ 


q11 -to  [ Thtz+Tdiode]  [T2 htz+T2 diode\  Ar 


(26) 
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Using  typical  experimental  values  for  Theater  and  Tdiode 
and  the  film  area  associated  with  the  heater  strip,  a 
resistance  value  for  radiative  heat  transfer  is  calculated 
between  the  heater  strip  and  the  cold  head  holder  using 


R 


radiative “ 


_  1  heater  1  diode 

q"‘A 


(27) 


and  the  resulting  value  of  Rradl-atjve  is  given  by 


R 


radiative " 


( I’htr  +  I'diode}  ^  ^  htr^^  diode ^  ^ 


(28) 


The  value  for  Rradiatjve  is  found  to  be  on  the  order  of 
104  K/W.  The  values  for  Rradiative  are  at  least  one  order  of 
magnitude  greater  than  the  values  for  Rsub<htr-  Consequently, 
radiation  losses  are  not  included  in  the  thermal  resistance 
network  or  error  analysis  considerations. 

5 .  Uncertainty  in  Rb 

The  values  of  Rb  are  plotted  in  Figure  22  against  Ts 
showing  error  bars  for  each  data  point.  As  anticipated  the 
error  in  each  value  is  minimal.  Having  used  cubic  spline 
interpolation  for  the  values  of  Th  and  Ts,  errors  due  to 
temperature  measurements  in  previous  studies  have  been 
eliminated . 
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c.  Rb  VARIATION  WITH  HEAT  FLUX 

Having  completed  the  error  analysis,  a  plot  of  Ts  vs 
the  corrected  values  of  Rb  is  shown  in  Figure  23  without 
error  bars.  Referring  to  the  plot,  the  values  for  Rb  with 
varying  heat  flux  are  shown  for  the  data  taken  at  average 
values  of  Ts  equal  to  91.2,  91.5,  110.8  and  113.6  K.  Values 
acquired  at  134  and  157.5  K  are  also  shown. 

It  is  clear  from  the  plot  that  the  majority  of  Rb 
values  over  the  entire  range  of  Ts  are  between  0.005  and 
0.001  W/cm2-K.  Despite  varying  the  values  of  heat  flux  at 
the  different  values  of  Ts,  no  discernable  dependency  of  Rb 
on  heat  flux  can  be  observed  from  this  plot.  Examining  Rb 
from  an  different  perspective,  heat  flux  vs  Rb  is  plotted  in 
Figure  24 . 

The  plot  of  heat  flux  vs  Rb  also  shows  no  overall  trend 
in  Rb  as  heat  flux  is  varied.  Given  that  the  variations  in 
heat  flux  in  this  study  were  done  to  explore  hysteresis 
effects,  it  is  apparent  that  the  individual  hysteresis  runs 
themselves  must  be  plotted  to  observe  any  trends  of  Rb  with 
heat  flux.  The  individual  hysteresis  runs  may  reveal  that 
just  as  important  as  the  current  value  of  heat  flux  is,  heat 
flux  history  of  the  sample  may  be  equally  important. 
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D.  HYSTERESIS  EFFECTS 

Cycling  of  the  heat  flux  at  fixed  values  of  Ts  was 
performed  to  determine  if  there  were  any  hysteresis  effects 
on  Rb .  Values  of  Ts  =  91.2,  91.5,  110.8,  113. 6K  were  used 
during  these  data  runs.  Plots  of  Heat  Flux  vs  Rb  were 
generated  and  are  shown  in  Figures  25  through  28.  Arrowheads 
are  used  on  the  plots  to  show  the  direction  in  which  heat 
flux  was  varied. 

Based  on  the  plots  alone,  the  conclusion  may  be  drawn 
that  Rb  does  display  some  sort  of  hysteresis  effect  when 
heat  flux  is  varied  to  a  maximum  and  back  to  a  minimum. 
However,  for  similar  values  of  Ts  (i.e.,  91.2  and  91.5  K  / 
110.8  and  113. 6K),  the  hysteresis  pattern  displayed  is  not 
consistent.  This  is  especially  noted  at  the  Ts  values  of 
91.2  and  91.5  K  where  the  values  used  for  heat  flux  were 
very  well  matched.  However,  with  all  the  hysteresis  runs, 
the  starting  value  of  Rb  is  always  greater  than  the 
finishing  values  of  Rb.  This  effect  is  believed  to  be  stress 
related  and  two  hypotheses  are  proposed. 

First,  it  has  been  shown  that  the  film  peeling  stress, 
that  is,  the  normal  stress  at  the  film/substrate  interface, 
becomes  more  negative  (i.e.,  more  compressive)  with 
increasing  q"  [Ref.  30] .  This  increase  in  compressive  stress 
causes  a  decrease  in  Rb,  as  seen  in  Figures  22  through  25. 

.As  q"  is  reduced  from  its  maximum  value,  some  of  that  high 
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compressive  stress  may  not  be  relieved,  leaving  the  film 
with  a  greater  compressive  (negative)  peeling  stress  than 
the  original  condition,  and  thus  a  lower  final  Rb  compared 
to  the  initial  Rb.  Second,  as  a  result  of  the  increasing 
compressive  stress,  microstructural  changes  are  produced 
near  the  film/substrate  interface.  These  are 
microstructural  changes,  such  as  the  movement  of  vacancies 
or  other  defects,  which  do  not  change  back  as  the  stress  is 
relieved.  These  changes  in  the  microstructure  will  affect 
the  microscale  heat  transfer  characteristics  between  the 
film  and  the  substrate  reflecting  in  the  changes  in  Rb 
observed . 

It  is  known  that  YBCO  displays  a  hysteresis  effect  on 
its  thermal  conductivity  upon  heating  and  cooling  [Ref.  30] . 
However,  the  presence  of  this  effect  cannot  be  determined 
from  the  hysteresis  runs.  What  is  clear  is  that  control  of 
variables  in  the  process  is  crucial  to  determining  the 
presence  of  hysteresis  effects. 

Current  experimental  techniques  used  in  this  study 
control  Ts  by  manually  cycling  the  cryo  cooler  on  and  off  in 
an  attempt  to  maintain  the  temperature  Tg.  Although  the  Lab 
View  program  will  not  acquire  data  until  steady  state  is 
reached,  duplication  of  conditions  from  one  data  run  to  the 
next  is  not  possible.  Considering  the  hysteresis  run  at  an 
average  value  of  Ts  =  91.2  K,  the  standard  deviation  of  Ts 
is  0.92  K.  This  results  in  a  variation  in  thermal 
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conductivity  of  YBCO  in  the  c-plane  from  2.146  to  2.138 
W/m-K  and  a  variation  in  the  thermal  conductivity  of  MgO 
from  277.27  to  290.86  W/m-K.  This  effect  resulting  from  poor 
control  of  Ts  is,  in  itself,  enough  to  mask  any  true 
hysteresis  effects  that  may  be  present.  This  serves  to 
emphasize  the  importance  of  controlling  the  experimental 
parameters . 

As  a  direct  result  of  this  study,  the  cryo  cooling 
apparatus  currently  being  designed  will  use  an  electric 
heater  to  control  the  temperature.  In  conjunction  with  the 
Lab  View  program,  the  resistance  value  of  the  sensor  will  be 
controlled.  A  constant  sensor  resistance  displayed  by  Lab 
View  represents  a  constant  value  of  Ts.  The  Lab  View  program 
additionally  allows  the  control  of  heater  current  and 
voltage  to  obtain  the  desired  power  input.  Thus,  having 
attained  positive  control  over  the  substrate  temperature  and 
the  heat  flux,  duplication  of  experimental  conditions  will 
be  possible.  Subsequent  data  runs,  no  matter  what  the 
purpose,  should  yield  more  useable  data. 
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E.  FINITE  DIFFERENCE  METHODS 

In  an  effort  to  predict  values  of  Rb  under  different 
experimental  conditions,  finite  difference  numerical 
techniques  were  used  to  model  the  physical  situation.  The 
control  volume  method  illustrated  in  Figure  29  was  used  in 
conjunction  with  a  graded  network  to  model  the  film  and 
substrate  [Ref.  31,  32].  Using  these  techniques,  two 
different  programs  were  written  using  Microsoft  Fortran 
Version  5.1  to  arrive  at  a  solution  using  the  Gauss-Siedel 
method:  one  based  solely  on  temperature  distribution  and  one 
using  both  temperature  constraints  and  heat  input  [Ref  33] . 

Focusing  only  on  the  temperature  distribution  within 
the  substrate,  a  graded  network  was  designed  as  illustrated 
in  Figure  30;  the  program  listing  is  contained  in  Appendix 
B.  Program  inputs  are  Th,  Ts  and  Tdiode  which  yield  outputs  of 
Rb  experimental,  Rb  calculated  and  temperature  distribution 
within  the  substrate.  The  grading  of  the 

network  was  controlled  by  programming  the  spacing  of  the 
nodes  in  both  the  I  and  J  directions.  Tolerance  between 
previous  node  values  and  calculated  values  is  also  inputted. 
The  program  initializes  all  nodes  to  Tdiode  and  then  imposes 
Th  and  Ts  on  the  network  and  subsequently  iterates  until 
convergence  within  the  prescribed  tolerance  between  the 
current  and  the  previously  calculated  values. 
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Figure  29.  Finite  Difference  Control  Volume 
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Figure  30.  Temperature  Distribution  Graded  Network 
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Figures  31  and  32  show  the  contour  plots  generated  for 
inputted  temperature  values  of  Th=92.8  K,  Tg=92.6  K  and 
T  =85.5  K  resulting  from  iteration  based  on  temperature. 

diode 

The  contour  plots  provide  a  qualitative  insight  into  the 
temperature  distribution  in  the  substrate.  Temperature 
contours  emanate  from  the  heater  strip  and  rapidly  drop  off 
to  values  approaching  Tdiode.  Consequently,  with  substrate 
dimensions  of  10  mm^  by  1  mm  thick,  the  heat  affected  zone 
is  very  small.  Figure  31  displays  a  substrate  thickness  of 
0.2  mm  from  the  heater  strip  and  Figure  32  displays  a 
substrate  thickness  of  0.04  mm  from  the  heater  strip. 

In  addition  to  the  temperature  distribution,  upon 
inputting  experimental  values  for  Th,  Ts  and  T^.^,  the 
temperature  difference  between  the  heater  strip  and  the 
substrate  directly  beneath  was  calculated.  Temperatures  for 
the  nodes  in  the  substrate  directly  beneath  the  heater  strip 
are  averaged.  The  temperature  difference  between  the  heater 
strip  and  the  substrate  was  calculated  using  Th  and  this 
average  value.  It  is  this  temperature  difference  that  the 
program  uses  with  the  heat  flux  to  calculate  the  values  of 
r  These  values  are  typically  within  1  order  of  magnitude 
of  the  actual  experimental  values,  but  do  not  accurately 
predict  the  experimental  values . 
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Vertical  position  (mm) 


Figure  31.  Substrate  Heat  Affected  Zone 
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Figure  32.  Temperature  Distribution  Close  to  the  Heater 


To  provide  a  mechanism  for  inputting  heater  power  into 
the  program,  the  temperature  distribution  program  was 
modified.  The  program  proceeded  as  before  iterating  on 
temperature  until  convergence  was  achieved.  Then  the 
constraint  of  holding  the  heater  temperature  at  Th  is 
removed  and  heat  generation  is  introduced  to  apply  the 
appropriate  power  to  the  heater  strip.  The  resulting 
adjustment  in  node  temperatures  yielded  values  of  Rb 
calculated  on  the  order  of  10'5  W/cm2-K,  which  is  far  less 
than  any  of  the  experimental  values  of  Rb  generated. 

It  is  believed  that  the  large  discontinuity  in  thermal 
conductivities  at  the  interface  of  the  heater  and  the 
substrate  results  in  temperatures  being  driven  to  lower 
values.  In  addition,  Microsoft  Fortran  running  on  a  PC  has  a 
limitation  of  approximately  60  by  60  nodes  of  memory.  This 
was  a  severe  limitation  because  the  program  results  were 
found  to  be  very  sensitive  to  grid  size.  The  memory 
limitation  precluded  performing  an  in  depth  grid  sized 
sensitivity  study .  In  an  attempt  to  overcome  the  averaging 
effect  on  the  thermal  conductivity  that  results  from  the 
finite  difference  method  of  solution  and  the  limitations  on 
grid  size,  a  new  program  was  written.  Given  that  the  heated 
affected  zone  is  small  compared  to  the  size  of  the  substrate 
and  that  the  sensor  strip  produces  no  joule  heating, 
symmetry  was  used.  A  centerline  was  passed  through  the 
heater  strip  and  only  half  the  substrate  was  modeled  as 
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illustrated  in  Figure  33.  The  program  functioned  in  the  same 
manner  as  the  previous  program,  however,  the  heater  strip 
was  modeled  by  numerous  layers  of  nodes  instead  of  just  one 
layer  of  nodes  in  an  attempt  to  overcome  the  averaging 
effects.  Despite  these  efforts,  the  new  program  still 
yielded  values  for  Rb  on  an  order  of  magnitude  of  10"5  W/cm2- 
K.  Once  again,  the  results  were  very  sensitive  to  grid  size. 

Poor  results  from  the  finite  difference  programs  are 
certainly  contributed  to  by  the  inability  to  use  a  small 
enough  grid  size.  However,  an  inability  to  model  factors 
such  as  microscale  heat  transfer  effects,  changes  in  peeling 
stress  of  the  thin  film  and  factors  included  in  the  acoustic 
mismatch  theory  probably  also  contributed  to  the  disparity 
between  experimental  and  calculated  values.  The  purpose  and 
the  premise  of  the  program  are  valid.  Consequently, 
restructuring  the  program  for  more  grid-sized  flexibility 
using  a  main  frame  computer  is  certainly  warranted. 
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V.  RECOMMENDATIONS 


This  study  has  focused  on  process  improvement.  The  use 
of  the  new  straight  line  etching  pattern,  elimination  of  the 
lock-in  amplifier,  use  of  a  pin  connector  and  the  Lab  View 
program  have  certainly  improved  the  data  acquisition 
process.  Utilization  of  the  cubic  spline  interpolation  for 
determining  heater  and  sensor  strip  temperatures  has 
eliminated  the  largest  source  of  error  from  previous  studies 
which  was  the  AT  calculation.  However,  further  improvement 
in  parameter  control  is  needed.  Utilization  of  the  heater  in 
the  cryo  cooler  to  control  the  resistance  of  the  sensor 
strip  and,  consequently,  substrate  temperature  is  essential. 

Once  substrate  temperature  control  is  achieved,  data 
runs  can  be  duplicated  to  a  reasonable  level  of  assurity. 

All  data  runs  from  this  study  should  then  be  redone. 
Particularly  the  hysteresis  runs  where  we  have  some 
indication  of  the  hysteresis  effects  being  present  but,  due 
to  poor  parameter  control,  the  results  are  inconclusive.  In 
addition,  it  is  recommended  that  fatigue  runs  be  done  as 
well.  Cycling  heat  flux  between  fixed  values  in  a  repetitive 
manner  may  indicated  whether  film  failure  or  progressive 
changes  in  Rb  will  occur  under  intermittent  service  use. 

It  is  recommended  that  the  finite  difference  program  be 
redone  on  a  main  frame.  This  would  remove  limitations  due  to 
computer  memory.  Additionally,  finite  element  methods  should 
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be  employed  if  possible.  Modeling  the  geometry  and  varying 
the  grid  size  could  be  done  with  much  greater  flexibility 
using  a  finite  element  program.  The  ability  to  model  the 
physical  situation  on  the  computer  is  important.  Given  the 
amount  of  time  in  the  lab  that  it  takes  to  acquire  data,  the 
program  could  be  a  valuable  tool  for  planning  data 
acquisition  runs. 

In  closing,  the  control  of  experimental  parameters  and 
computer  modeling  should  be  used  in  the  performance  of 
further  heat  flux  variation  studies.  Re-running  past 
experiments  must  be  done  to  establish  a  consistent  record  of 
Rb  variation.  In  addition,  further  areas  of  heat  flux 
variation,  such  as  fatigue  affects  must  be  studied. 
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APPENDIX  A.  SUBSTRATE  TEMPERATURE  DISTRIBUTION 


Referencing  Figure  14  and  using  the  heat  conduction 
equation  for  two  dimensional  steady  state  heat  transfer  we 
get  [Ref .  34] , 


d2  t  +  a2  T 
dx 2  dy2 


(1) 


with  the  boundary  conditions  (BC) : 


(i) 


@  x=0 


dx 


(ii) 


@x=a 


dT 

dx 


=0 


(iii) 


@x=a  T= 0 


(iv) 


@y=b 


dT  X 
—  =f{x) 

dy 


where  f (x) = 


0<x<x, 
k  1 

o  xx<x<a 
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Boundary  condition  (iii)  is  defined  for  Figure  14  to  make 
the  problem  homogeneous : 

T=T-Tdiode  (2) 

where  T  is  the  actual  temperature  distribution.  Since  we 
are  interested  in  the  actual  temperature  distribution  (AT) , 
equation  2  can  be  ignored  later. 

The  solution  for  equation  (1)  is  given  by 

T=X-Y 


where 

X=C1cos3x+C2sin(3x  (3) 


Y=C3cosh(5y+C4sinh|3y  (4) 


such  that 


—  =-C1 3sin(3x+C2(3cos3x 
dx 


(5) 


- =-C,(3sinh3y+C.cosh(3y 

dy  J 


BC  (i)  and  (5)  gives 


=  0 


(6) 
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BC  (ii)  and  (5)  gives 


-C1PsinPa=0 


=> 


P n=I™  =>  Pn  = 


ni i 
a 


From  BC  (iii)  and  (4) 


C3=0 


(7) 


(8) 


(9) 


Using  the  orthogonality  property  we  get 

faf(x)  cosp„xdx=  f  *CnPncoshPnb cos2Pnxdx 
Jo  Jo 


LHS=  f*1-? COS$nxdx+  f  a0  -cos  P^xdx 


Hence 


T=Cncos  pnxsinhpny 


and  BC  (iv)  and  (6)  give 


H =caPacos  P„X  Sinhpny=f  (x) 
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i?ffi?=|cnpnc  OShp^  (11) 

Equating  the  LHS  to  the  RHS 

■^-sinpj2xl--|capjlcoshpa7 
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Thus 


,  2gsinPnx1 
icp„acoshPny 


and 


r=£  c„cospnxsinhPny 

D= 1 


T= 


2gsinPJ:x1  .cospnXS^n]1p^y 
k$2nacosh$ny 


The  temperature  difference  between  the  area  below  the 
strip  and  the  sensor  strip  is  given  by  the  difference 
between  the  mean  temperatures  under  each  strip. 

The  mean  temperature  below  the  heater  strip  is 


xi 


L Imean 


m 

-Ei 

10  =  1 


2gsinpj,x1sinp„x] 


L imean 


=E 

n=l 


2gsinpnx1sinpnx1^ 
k^2nax1  0 


(12) 


(13) 


heater 


(14) 
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The  mean  temperature  below  the  sensor  strip  is 


T 

^  2  mean 


1 

(a-x2) 


[ 


2gsinpnx1sinp/Ix  * 

^PnaPn 


Llmean 


=-E 


2qsin^nx1sin^nx2 


n=i  k$na  (a-x2) 


(15) 


The  temperature  drop  is  then  calculated  as 

A  T=  T  —  T 

Imean  **  2  mean 


A T=y  2gsinpnx1sinpnx1  A  2gsinPflx1sinpnx2 
n=i  JcPfl&jq  n=i  .kp^a  (a--x-2) 


Ar=f'  2gsinB”xi  (  sinPi?xi  +  sinPn*2 ) 
hi  £p^a  Xx  a-x3 
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The  series  converges  and  is  approximated  for  n=3 .  From 
equation  (7)  and  the  values  of  the  constants  from  Figure  14 
equation  (16)  becomes 


57=0 . 0154  — 
k 


(17) 


Using  the  experimental  values  from  each  data  run  in 
conjunction  with  the  value  of  Ts  a  corrected  value  of  Rb  is 
calculated . 
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APPENDIX  B.  FINITE  DIFFERENCE  PROGRAM 

TEMPERATURE  DISTRIBUTION  OF  A  THIN  FILM  ON  A  SUBSTRATE 
SOLUTION  BY  GAUSS -SEIDEL  ITERATION  METHOD 

C  INTEGER  I , J, LI , LJ 

INTEGER  L 

PARAMETER ( LI = 5 7 , L J= 5 9 ) 

DIMENSION  PX (LI, LJ) ,  PY(LI,LJ),  T(LI,LJ)  , TTEMP (LI , LJ) 
REAL  DISTX,  DX,  DISTY,  D2 ,  D3 ,  XA,  YA,  DIFF,  TOLL, 

KS ,  KF 

OPEN (UNIT=10 ,  FILE= ' CONTOUR ' ,  STATUS= ' UNKNOWN ' ) 

C  DEFINE  TOLERANCE 

TOLL=0 . 005 

P=  .182383 

Th=92 . 837 

Ts=92 .672 

Tcf =85 . 55 

X= (Th+Ts) / 2 

KF=2 .605-0. 00506 *X 

KS=22 83 .61772-39. 9669583*X+0 . 245369913 *X**2 
-0 . 0 005 08 78108 *X* *3 

C  DEFINE  X  AND  Y  NODE  POSITIONS 

C  DEFINE  1ST  ROW  X  POSITIONS 

DISTX=0 .0 
DO  10  J=1 , LJ 

C  DEFINE  TOP  ROW  SPACING  USING  IF  STATEMENTS 

IF ( J . LE . 15 )  DX=0 . 01 

IF( (J.GE.16) .AND. (J.LE.49) )  DX=0.001 
IF ( J . GT .49)  DX=0 . 01 


C  WRITE  PX  AND  PY  FOR  J  POSITION 

PX ( 1 , J ) =DISTX 
PY ( 1 , J ) =0.0 

C  ESTABLISH  SPACING  FOR  THE  J+l  POINT 

DISTX=DISTX+DX 
10  CONTINUE 

C  EXPAND  1ST  ROW  POSITIONS  INTO  Y  POSITIONS 

DO  20  1=2, LI 
DO  30  J=l. LJ 

PX ( I , J ) =PX (1 , J) 

IF ( I . EQ . 2 )  DISTY=0 . 00001 
IF ( I . EQ . 3 )  DISTY=0. 00001 
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on  o  nm  noon 


IF ( I . EQ . 4 )  DISTY=0 .00001 
IF ( I . EQ . 5 )  DISTY=0 . 00001 
IF ( I . GT . 5 )  DISTY=0 . 01 
IF ( I . GT . 51 ) DISTY=0 . 1 
PY ( I , J) =PY ( I - 1 , J) +DISTY 
30  CONTINUE 

20  CONTINUE 

C  INITIALIZE  T ( I , J)  THE  TEMPERATURE  ARRAY 

DO  50  1=1, LI 
DO  60  J=1 , LJ 
T ( I , J) =Tcf 
60  CONTINUE 

50  CONTINUE 

DO  52  J=16 , 31 
T (1, J) =Th 

52  CONTINUE 

DO  53  J=36 , 49 
T (1, J) =Ts 

53  CONTINUE 

USE  GAUSS -SEIDEL  ITERATIVE  TECHNIQUE  TO  CALCULATE  THE 
TEMPERATURES 

SET  THE  ITERATION  COUNTER  N=1  FOR  THE  FIRST  ITERATION; 
SET  THE 

.  TOLERANCE  COUNTER  K=0 
N=1 
L=0 

PRINT*, ' L= '  ,  L 
K=0 

CALCULATE  VALUES  FOR  THE  HEATER  STRIP 
IF  (L.EQ.0)  GOTO  59 

DO  57  J=16 , 31 
D3=PY (2 , J) -PY (1, J) 

XA= (PX(1, J+l) -PX(1, J-l) ) /2 

compensate  power  for  the  size  of  the  control  volume 
q=P*1000* (0.1/. 16) 

KS=2893 . 84258-657. 933* (92 . 754-30) **0.337 
KS=100*7762 .715*92.754** (-1. 6941479) 

KTOT= (KS+KF) / ( KS  *  KF ) 

TTEMP (1, J) =T (2 , J) + (q*D3*KTOT) / (XA*2) 

DIFF=ABS (TTEMP (1 , J) -T ( 1 , J) ) 

IF( (DIFF-TOLL) .GT.0.)  K=K+1 
T (1 ,  J) =TTEMP (1 , J) 
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57  CONTINUE 

C  UPPER  LEFT  HAND  CORNER  NODE 

59  PRINT*, 'GO  TO  59' 

D2=PX (2,2) -PX (2,1) 

D3=PY (3,1) -PY (2,1) 

XA=D2/2 . 

YA=D3 / 2 . 

TTEMP (2,1)= (YA*D3*T (2,2) +XA*D2*T (3,1))/ (YA*D3+XA*D2) 
DIFF=ABS (TTEMP (2,1) -T(2,l) ) 

IF( (DIFF-TOLL) .GT.O. )  K=K+1 

T ( 2 , 1 ) =TTEMP (2,1) 

C  TOP  EDGE  NODE 

DO  70  J=2 , LJ-1 

IF(  (J.GE.16)  .AND.  (J.LE.31)  )  GOTO  75 
IF( (J.GE.36) .AND. (J.LE.49) .AND. (L.EQ.O) )  THEN 
GOTO  75 
END  IF 

D2=PX (2 , J+l) -PX (2 , J) 

D3=PY (3 , J) -PY (2 , J) 

D4=PX (2 , J) -PX (2 , J-l) 

XA= (D2+D4 ) / 2 . 

YA=D3 / 2 . 


TTEMP (2, J) = (YA*D3*D4*T (2 , J+l) +XA*D2*D4*T (3 , J) 
+YA*D2*D3*T (2 , J-l)  ) 

+/ (YA*D3  *D4+XA*D2  *D4+YA*D2  *D3 ) 

GOTO  65 

75  1=2 

D1=PY (I , J) - P Y ( I - 1 , J ) 

D2=PX(I, J+l) -PX ( I , J) 

D3=PY (1+1 , J) - PY ( I , J ) 

D4=PX (I , J) -PX(I, J-l) 

XA= (D2+D4 ) / 2 
YA= (D1+D3 ) / 2 

C  KF=5 . 0 

C  KS=2893 . 84258-657 . 933* (92 . 754-30) **0.337 
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TTEMP ( I , J) = (2*XA*KF*D2*D3*D4*T(I-1, J) 
++YA*D1* (KS+KF) *D3*D4*T(I,J+1) 

++XA*D1* (KS+KF) *D2*D4*T(I+1, J) 

++YA*D1* (KS+KF) *D2*D3*T(I, J-l) ) / 

+  ( 2  *XA*KF*D2  *D3  *D4+YA*D1* (KS+KF) *D3*D4 
++XA*D1* (KS+KF) *D2*D4+YA*D1* (KS+KF) *D2*D3) 

65  DIFF=ABS (TTEMP (2, J) -T(2,J) ) 

IF( (DIFF-TOLL) . GT . 0 . )  K=K+1 

T  (2  ,  J)  =TTEMP  (2  ,  J) 

70  CONTINUE 

C  UPPER  RIGHT  HAND  CORNER  NODE 

D3=PY  (3  ,  LJ)  - PY  ( 2  ,  L J ) 

D4=PX  (2  ,  J)  -PX(2,  J-l) 

XA=D4/2 . 

YA=D3 / 2 . 


TTEMP (2, LJ)  =  (XA*D4  *T  (3 , LJ) +YA*D3 *T ( 2 , LJ-1 ) ) / (XA*D4+YA*D3 ) 
DIFF=ABS (TTEMP (2, LJ) -T(2,LJ)  ) 

IF ( (DIFF-TOLL) .GT. 0 . )  K=K+1 

T ( 2 , L J ) =TTEMP ( 2 , L J ) 

C  SET  UP  DO  LOOP  TO  CALCULATE  LH,  INTERIOR  AND  RH  NODES 

FROM 

C  1=3  (ROW  3)  UNTIL  I=LI-1  (NEXT  TO  LAST  ROW) 

DO  80  1=3, LI-1 

C  LEFT  HAND  SIDE  NODES 

D1=PY (1,1) - PY (1-1,1) 

D2=PX (1,2) -PX(I, 1) 

D3=PY (1+1, 1) - PY (1,1) 


XA=D2/2. 

YA= (D1+D3) /2 . 


TTEMP (I, 1)  =  (XA*D2*D3*T (1-1, 1) +YA*D1*D3  *T ( 1 , 2 ) 
+XA*D1*D2  *T (1  +  1,1) ) / 

+ (XA*D2*D3+YA*D1*D3+XA*D1*D2) 

DIFF=ABS (TTEMP (I, 1) -T (1,1) ) 

IF ( (DIFF-TOLL) .GT.O.)  K=K+1 
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T (I , 1) =TTEMP (I, 1) 

C  INTERIOR  NODES 

DO  90  J=2 , LJ- 1 

D1=PY ( I , J) -PY (1-1 , J) 

D2=PX ( I , J+l) -PX ( I , J) 

D3=PY ( 1+1 , J) - PY ( I , J ) 

D4=PX ( I , J) -PX(I, J-l) 

XA= (D2+D4) /2 . 

YA= (D1+D3 ) / 2 . 

TTEMP (I, J)  =  (XA*D2*D3*D4*T (1-1, J) +YA*D1*D3  *D4  *T ( I , J+l) 
++XA*D1*D2*D4*T ( 1+1 , J) +YA*D1*D2*D3  *T ( I , J-l) ) / 

+ (XA*D2*D3*D4+YA*D1*D3*D4+XA*D1*D2*D4+YA*D1*D2*D3 ) 

DIFF=ABS (TTEMP (I , J) -T ( I ,  J)  ) 

IF ( (DIFF-TOLL) .GT. 0 . )  K=K+1 

T ( I , J) =TTEMP ( I , J) 

90  CONTINUE 

C  RIGHT  HAND  SIDE  NODES 

D1=PY ( I , LJ) -PY (1-1,  LJ) 

D3=PY ( 1+1 , LJ) -PY ( I , LJ) 

D4=PX ( I , LJ) -PX ( I , LJ-1) 

XA=D4/2 . 

YA= (D1+D3 ) / 2 . 

TTEMP (I,LJ) = (XA*D3*D4*T (I-1,LJ) +XA*D1*D4*T (I+1,LJ) 
++YA*D1*D3  *T ( I , LJ-1) ) / (XA*D3*D4+XA*D1*D4+YA*D1*D3) 
DIFF=ABS (TTEMP ( I , LJ) -T ( I , LJ) ) 

IF  (  (DIFF-TOLL)  .GT. 0 . )  K=K+1 

T ( I , L J ) =TTEMP ( I , L J ) 

80  CONTINUE 

IF  ( (K.EQ.O) .AND. (L.EQ.l) )  THEN 
PRINT*, ' T (2 , 25) =' , T (2 , 25)  , '  K='  ,K 

GO  TO  100 

END  IF 

IF (K.EQ.O)  THEN 

L=1 

N=1 
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PRINT* GOTO  55' 

C  GO  TO  55 

GO  TO  100 

ELSE 

PRINT*, ' T ( 2 , 25 ) = ' , T (2 , 25 )  ,  '  K='  ,K 

N=N+1 
GO  TO  55 
END  IF 


100  DO  105  1=1,26 
DO  110  J=1 , LJ 

WRITE (10,102) I , J, PX (I , J) , 1-PY ( I , J) ,T(I, J) 
102  FORMAT ( IX, 14,14, 3F12 . 8 ) 

110  CONTINUE 
105  CONTINUE 

PRINT*, ' L= ' , L 
PRINT*, ' N= ' , N 

Th=0 

DO  107  J=16 , 31 
TAVG=Th+T ( 1 , J) 

Th=TAVG 

107  CONTINUE 
Th=TAVG/l6 
PRINT*, ' Th= ' , Th 

Tc=0 

DO  108  J=16 , 31 
TAVG=Tc+T (2 , J) 

Tc=TAVG 

108  CONTINUE 
Tc=TAVG/l6 
PRINT*, ' Tc= ' , Tc 

Rb= (Th-Tc) / (P/0.0016) 

DELT=Th-Tc 

PRINT* , ' DELTA  T= ' , DELT 
PRINT* , ' Rb= ' , Rb 
Rbexp= (Th-Ts) / (P/0 . 0016) 

PRINT* , ' Rbexp= ' , Rbexp 

PRINT*, ' KS= ' , KS 

END 
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